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SIMIARY 


Linear theory is used to caJLcvilate the surface pressures, field 
pressvires, and wave dreig of nonlifting humped and Indented hodles of 
revolution in supersonic flow end the results are compared with the cor- 
responding properties of a smooth basic body. The calculations show 
that relatively small surface irregularities cause large pressure dis- 
turbances both on the body and in the field. Application of a correction 
for the curvature of the characteristics substantially changes the nature 
of the pressure distribution, particularly in the disturbance field gen- 
erated by a stirface irregularity. The drag results indicate that, in 
general, the wave drag increases if volume is added to or subtracted from 
a smooth basic body so as to produce a surface irregularity on the basic 
body. The total wave drag of nonsmooth bodies consists of two relatively 
simple terms - the drag of the basic body and the drag of the body deter- 
mined from the area distribution of the bump or indentation alone - and 
a complicated interference term. The Interference parameters are pre- 
sented for two body types and a range of valiaes is suggested for use in 
obtaining a reasonable estimate of the wave drag of nonsmooth bodies. 


INTRODUCTION 


The aerodynamic characteristics of bodies of revolution at super- 
sonic speeds have been considered in numerous investigations and several 
methods have been developed to describe the flow past such shapes. Tbs 
method of characteristics and the linear theory approximations of 
Von Karman (ref. l), Llghthill (ref. 2), and Ward (ref. 5) have proven 
especially usefvil in this regard. Although many calculations have been 
made by using these methods, most have been limited to smooth bodies. 

With the advent of the area rule, however, there has been considerable 
interest in nonsmooth shapes. 

The purpose of the present paper is to consider the aerodynamic 
properties of a simple class of nonsmooth bodies - namely, bodies of 
revolution which have irregularities in area distribution. Linear theory 
is used to calculate the surface pressures, field pressinres, eind wave 
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drag for nonliftlng bumped, and. indented, parabolic bodies of revolution 
in supersonic flow and the results are ccmpared with the corresponding 
properties of a smooth basic body. The Whitham theory (ref. 4) is used 
to correct the field, pressures for the curvature of the characteristics. 


SYMBOLS 


A 

B 

Cp 

D 


cross-sectional area 

interference -drag parameter appearing in equation (l4) 

p - p 

pressure coefficient, — ^ 


wave drag 


F(y) 



l/y - I 


f source-distribution function 


tr _ 7 + 1 


p3/2 


I 

M 

N 

P 

Fc 

<1 

R 

U 

V 


length of soTurce distribxrtion 
free-stream Mach nimiber 
soTirce-distribution constant 
local pressure 
free-stream pressure 


dynamic pressure. 



body radivis 
free-stream velocity 
volimie 
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axial eind radial disfances, respectively, in cylindrical 
coordinates 

coordinates defined by equations (l2) 
origin of bump or indentation source distribution 

characteristic number, x - pr 
3 = - 1 

A incremental value 

7 ratio of specific heats 

I dimimy variablje of integration 

p free-stream density 

0 disturbance velocity potential 
Subscripts : 

1 basic body 

2 bump or indentation 

SH Sears -Eaack body 

X partial derivative with respect to x 

r partial derivative with respect to r 

Primes denote derivatives with respect to the indicated argument 
of the function. 


x,r 

5,r 

y 


METHOD OF CALCULATION 

The linearized differential equation for axlsymmetrlc supersonic 
flow is 


P^xx - 0rr - f 9^ = 0 


( 1 ) 
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where ^ M Is the stream Mach ntnnher, and (j> Is the dls- 

turhance velocity potential. (The coordinate system Is defined In 
fig. 1.) A solution of equation (l) for flow past a body of revolution 
can he expressed as 



\/(x - 1)^ - 


( 2 ) 


where f(|) represents a dlstrlhution of soiorces along the axis of the 
body. The source distribution f(|) is related to the body configura- 
tion by the boundary condition of zero normal velocity at the body sior- 
face. With the aid of mass-flow considerations (see ref. 5)^ this 
boundary condition can be written as 


itR^(x) 




x-pR(x) 


(x - |)f(|)d| 

)/(x - - 3%^(x) 


(3) 


where R(x) is the body radius at any station x. 


Inasmuch as the inverse of the integral equation (eq.. 5 ) for the 
sovirce strength is not known, the method of analysis adopted was to 
consider the sovirce distribution as the primary variable. With this 
approach, nonsmooth bodies of revolution can be generated by the super- 
position of two known source distributions. The source distribution 
fl(x) gives rise to the basic body which lies between x - 0 and 

X = Z]_. The bump or indentation is created by superposition of the 
source distribution f 2 (x) on fQ^(x) from x = Xq to x = Xq + Ig, 
where Xq denotes the origin of fg and Zg is the length of fg. 
Thus, 


f (x) = fi(x) 

(0 ^ X ^ Xq j 

f (x) = fi(x) ± f 2 (x) 

+ 

VII 

VII 

f (x) = fi(x) 

fxQ + Z2 < X S 


(4) 


A plus sign before corresponds to a bump and a min\is sign, to 

an indentation. 
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The expressions for f^Cx) and f2(x) used in the analysis are 


fl(x) = + 2x^) 


and 




=22 


^2^(x - Xq) - 3 ^ 2 (^ - + 2(x - Xq)' 


( 5 ) 


The functions f^^Cx) and f2(x) each give rise to a parabolic body of 
revolution in the slender-body approximation 


f (x) = A'(x) 


( 6 ) 


where A(x) is the cross-sectional area distribution and the prime 
denotes the derivative with respect to the argument of the function. 

The parameter N^_ or N2 which appears in equations (5) can be inter- 
preted geometrically as the slender-body approximation (eq. 6 ) to the 
fineness ratio of the body generated by f]_ or f2 alone. 


Calculation of Body Shape 

Equation (3) may be written in terms of the source distributions fj_ 
and f2 as 


itR^(x) 


-fo 


x-^R(x) (x - 


(x - 


X 


x-pE(x) (x - Ofp(Od| 


(x - 1 )^ - p%^(x) 


(7) 

The calculation of body shape from this equation is tedious Inasmuch as 
R(x) occurs in each term. Consequently, in the foUoirLng paragraphs 
the equation is simplified and expressed in a form which partleilly sepa- 
rates the calculation of bump shape frcm the calculation of basic body 
shape. 
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For certain types of problems, it is well known that "body shapes 
calculated hy equation ( 3 ) can he closely approximated, hy preassigning 
a value to R(x) on the right-hand side of the equation. If R(x)/z 
is pmial 1 compared with unity, R(x) may he set equal to zero and the 
equation reduces to the simple slender-hody relation (eq. 6) . For 
quasl-cyllndrlced. configurations (ref. 2), R(x) may he assigned a 
constant value on the right-hand side of equation (3)» In reference 6, 
which considers the problem of hoattall bodies of revolution having 
TTiiniTmiTn drag, R(x) is assigned the values corresponding to the cone 
that passes through the nose and base of the body. 

,^-dIhe approximate methods outlined suggest that some reasonable sim- 
plifications can be made in equation ( 7 )- For the body shapes considered 
herein, R(x) is veiy smal 1 compared with the length of fj_ and, con- 
sequently, the slender-body approximation can be made in the first inte- 
gral of equation ( 7 )- However, inasmuch as R(x) is not necessarily 
small compared with the length of f£, some other approximation is 

required in the second integral. An approximation that has the merit 
of simplicity £ind accuracy is to replace R(x) in the second Integral 
of equation ( 7 ) with the radlois of the basic body as determined from 
the slender-body relation 2 jiRjl(x)R 3_' (x) = f^^Cx). Then 


rtR^(x) = f fi(|)d| ± f 

do d y 


x-pR^(x) (x - 




( 8 ) 


where the first integral, represents the area distribution of the slender 
basic body and the second Integral, which for a given Mach number depends 
only on R^^Cx) and f£, represents the area distribution of the surface 

irregularity . 

A comparison of the bimip area distributions AA. as obtained from 
equations (6), ( 7)7 a nd (8) for = 12, N 2 = 8, Xq^Zj^ = O. 5 O, 

Z 2 /Z 1 = 0 . 25 , and M = 1 / 2 " is presented in figure 2 and a comparison 

of the total area distribution is shown in figure 3- approximation 

afforded by equation (8) for both the bump area and the total area is 
considered satisfactory, whereas the slender-body approximation (eq. 6) 
is considered Inadequate for the calcxolation of nonsmooth shapes. 


Calculation of Pressure Field 

T.inear theory .- The .linearized pressure coefficient at any point 
in the field is given by 


Cp - -20X - 0r^ 


(9) 
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where 


and 


!<:c 


1 . f'({)dt 

^^0 ]j{x - - P^r2 


[ = _L. (x - |)f '(|)d| 

^"'O |/(x - 1)2 - p2r2 


( 10 ) 


( 11 ) 


Whlthani*s correction .- It is well known that linear theory does 
not provide a first-order approximation to. the flow field at supersonic 
speeds that the approximation hecomes worse as the distance from 
the body Is Increased. This fact led Whltham (ref. 4) to develop a 
method for correcting the linear flow fields so that a first approxima- 
tion Is obtained everywhere. Whitham’s basic hypothesis Is that linear 
theory gives the correct valties for the velocity components but locates 
them Inconectly. By taking the curvature of the characteristics Into 
accoimt, Whltham deduces that the values of ^^(x,r) and 0j,(x,r) 

should be located at the coordinates (x,f ) where 


r = r 

X = pr - KB’(y)Vr' + y 


(12) 


In equation (l2) 


y = X - 3r 



f'(l)d| 


and 


K = y + ^ 

p5/2 
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When the flow-field data are relocated by means of equation {l2), 
it is found that, in some cases, the characteristic lines cross. Whitham 
interprets such intersections as indications of the presence of shock 
waves. Reference 4 presents further details and an interesting graphical 
method of locating the shock waves from the F(y) curve. 


Calculation of Wave Drag 

The wave drag of a distribution of sources along an aois can be 
related to the mcmientum transfer through a cylindrical surface enclosing 
the axis. The ejjpression for the wave drag is then obtained (see ref. 1, 
for example) as 


23tD 

q. 


~ Iq ^ Iq ^ f'(x)f'(Ologelx - dx d| 


(13) 


where D is the wave drag and q is the dynamic pressvire . When the 
source distributions given by equation (4) are substituted into equa- 
tion (13) an^ indicated integrations performed, the wave drag can 
be written as 


or 






+ 



^1 


(14) 


where D]_ is the wave drag of the basic source distribution fj_ alone, 
I >2 is the wave drag of f 2 alone, and D is the total wave drag which 
results from the combination of f^^ and f 2 » A plus sign preceding the 

second term on tne right-hand side of equation (l4) corresponds to a 
bump on the basic body and a minus sign corresponds to an indentation. 
The factor B appearing in the equation is a ccanplicated function 
of Xq/Z]_ and ^2/^l‘ 
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EESUIZES AND DISCUSSION 
Body Shape 


Scale drawings of the basic, humped, and. indented, bodies considered 
herein are presented in figure These shapes were calculated from the 
basic mass-flow boundary condition (eq.. ( 7 )) for Nq_ = 12 , N 2 = 8 , 

= 0 . 50 , ^ 2/^1 “ sj^d M = The differences between the 

basic body shape and the shapes of the humped and indented bodies appear 
finial 1 for the case selected. In order to illustrate more clearly the 
magnitude of the area change introduced by the surface irregularities, 
the cross-sectional area distrihirtions of these bodies are presented in 
figure 5 * 

The variation of hump area distribution with the strength of the 
hump source distribution is shown in figure 6 , where bump area distri- 
butions are presented for values of N 2 = 6, 8 , and 10. The Indentations 

for corresponding values of N 2 are slightly different frcm the bump 

area distributions shown Inasmuch as body shapes calculated from equa- 
tion ( 7 ) are inf licenced by the differences between the radii of the 
bumped and indented bodies. If equation ( 8 ) is used to calciiLate body 
shape, however, the area distributions of corresponding buTnpp nnfl inden- 
tations are the same. 


Pressure Field 

Ole siarface pressure distributions for the basic, bumped, and 
indented bodies corresponding to values of N-j_ = 12, = O. 5 O, 

^ 2/^1 = M = \f 2 ^ and N 2 = 6 , 8 , and 10 are presented in fig- 

ure 7 . Although the surface irregularities considered make only emp.n 
visual changes in body shape (for example, see fig. X) , the surface 
pressure disturbances caused by the irregularities are seen to be rela- 
tively large. It should be pointed out, however, that viscous effects 
would probably reduce the magnitudes of the pressure distiorbances arising 
from such smn. 1 l surface irregularities . 

In the limit as the bimip length goes to zero, the change in surface 
pressure coefficient due to the bump can be obtained exactly from the 

two-dimensional pressure relation ££rg = ^ AR', where AR' is the change 
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in body slope due to f 2 - A eomparison of tbe -valijes of at tbe 

surface as obtained from equation ( 9 ") and the two-dimensional relation 
is presented in figure 8 for the bump considered in the present paper. 

Field-pressure coefficients, with, and without the Whitham correc- 
tion, for the basic, bumped, and indented bodies (fig. 4), are presented 
in figures 9 } 10, arid U, respectively. Pressure distributions for each 
body shape are given at the four field locations = 0 . 125 , 0 . 250 , 

0 . 500 , and 1.0. 

The shift of the field pressures predicted by the Whdtham theoiy 
for the basic body (fig. 9 )) althou^ significant, is not large at the 
field locations shown, and the predicted shock at the nose of the body 
affects only a Rmall portion of the field. However, the VJhltham cor- 
rection in the region of the bump or Indentation (figs. 10 and U) is 
relatively large at an field stations represented, and the predicted 
shocks have a marked effect on the pressure distribution. 

In order to illTistrate the natinre of the calculated flow field about 
the basic, bimiped, and indented bodies, sketches of the Whitham corrected 
characteristics and shock locations are presented in figure 12. In the 
region of the surface Irregularity the corrected characteristics exhibit 
considerable ciirvature, whereas the 1 1 ne ar characteristics would have no 
cuirvature ant^ would point downstream at an angle of ^5° from the boi3y 
center line. 


Wave Drag 


The variation of D^D]_ with as calcvilated from equation (l4) 

for both bumped and indented bodies is presented in figure I 5 for 
^ 2/^1 “ 0*25 with Xq /^1 = 0 and 0.375^ and for ^ 2/^1 “ 0*50 with 
Xq|Zi = 0 and 0 . 25 . The values of chosen bracket the range 

of D^D]_ for the given value of for other 

values of and the given ^ 2/^1 between the curves shown. 

Some of the curves correspond to two locations of the source distri- 
bution f 2 : For. example, the bumped-body curve for ^ 2 . 1 ^ 1 . = 0*25 and 

Xq/Zi = 0.75 is Identical to that for Xq^Z-^ = 0. It should be noted, 

however, that the drag of a given nonsmooth body is not exactly the 
same in a forward and reverse flow because the bump area depends upon 
both the b\mip source strength and the configuration of the basic body. 

(See eq.. ( 8 ).) 


Examination of figure 13 shows that in most cases the wave drag 
Increases if volume is either added to or subtracted from a smooth basic 
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body 80 as to prod.uce a surface irregularity on the basic body. For 
example, the bumped, and, indented, bodies (fig. 4 ) have 4 l percent and 
22 percent more drag, respectively, than the smooth basic body. 

According to figure I3,, the ratios %/^2 

mary importance in the determination of the drag, whereas the location 
of the surface irregularity is of secondary in^jortance. The latter 
conclusion might have been anticipated frcm a result of Robert T. Jones 
presented in reference 7 that, for a Sears-Haack basic body, the drag 
Increment duj^ to adding or subtracting a bump is independent of the 
location of the bump. If fgg denotes the source distribution which 

would give rise to a Sears-RjEiack body in the slender-body approximation 
and f2 denotes any other source distribution shorter than fgg, Jones 

found that the drag of the shape given by superposition of f g on fgg- 
can be expressed as 


D 

®SH 


1+2 




(15) 


In equation (15), Dg denotes the drag of fg alone and denotes 

the drag of fgg- alone. The symbols Vg a-^d Vgg- are to be inter- 
preted in a manner similar to that for N for the parabolic shapes; 
that is, Vg (or Vgg^ is the voltime of the body defined by equation ( 6 ) 

for fg (or alone. 

The expressions for the drag given. by equations (l 4 ) and (15) nat- 
urally separate into three parts - the drag of the basic body alone, 
the drag of a body determined by the volume of the bimip or Indentation 
alone, and an Interference drag. In general, the drag of the separate 
components can be calculated quite easily, whereas the computation of 
the interference drag is rather involved. It is noteworthy that Jones 
discovered a case for which the interference drag is given by such a 
simple expression. 

^Jhen equation (15) is written in the form of equation (llf) for two 
Sears-Haack source distributions, the interference-drag parameter B 
is given by 
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These -values of B are shown in figure -l4 for comparison with 
those appearing in equation (l4) for the parahollc shapes considered in 
this paper. There is general agreement be-tween the magnitudes of B 
for -fche -various length ratios. This agreement suggests the possihillty 
of oh-baining a rapid estima-te of the dorag of nonsmdoth shapes hy selec- 
tion of a -value of B from figure l^J- on the basis of length ratio alone. 


CONCLUDING KE34ARKS 


Thp. supersonic pressure field and wave drag of nonlifting bumped 
and Indented bodies of revolixbion are calcula-fced by linear theory and 
the results are compared with the corresponding properties of -a smooth 
body. A compeirison of the calculated pressure distributions indicates 
that relatl-vely sma ll surface iiregularities cause large pressrire dis- 
-turbances both on the body and in the field. Application of a correc- 
tion for the cu 2 Tvature of tlie charac-fceristlcs subs-bantially changes the 
na-ture of the pressinre distribution, particularly in the disturbance 
field genera-fced by a surface irregularl-ty. Th e drag results Indicate 
that, in general, the wave, drag increases if -volume is added to or sub- 
tracted from a smooth basic body so as to produce a surface irregularl-fcy 
on the basic body. The total wa-ve drag of nonsmooth bodies consists of 
•fcwo relati-vely simple -terms - the drag of the basic body and the drag 
of the body de-termlned from the cross-sectional area distribution of 
the bump or inden-tation eilone - and a complica-ted interference term. 

The interference parameters are presented for -two body types and a range 
of -values is suggested for use in obtaining a reasonable estimate of the 
drag of nonsmooth bodies. 


Langley AeronautlcsLl Laboratory, 

National Ad-vlsory Commi-btee for Aeronautics, 
Langley Field, Va., May 1^, 1956. 
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Figure 2.- Ccai^jarlson of c han ge a in hiocly area introduced "by a bump source 
distribution as calculated from equations (6), ( 7 )^ ^d (8). = 12^ 
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= O.pO; ^ 0.25; M = \/2. 
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(a) Baeic and tiuraped 'bodies. 
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(b) ^ = 0.250. 

Figure 9‘- Continued. 
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(a) Basic tody. 

Figure 12.- Characteristic field and shock locations predicted hy Whitham 
theory for hasic^ buinped^ and indented "bodies of figure 4. Mm sfz. 
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(b) Bundled "bocly. 
Figure 12.- Continued, 
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(a) Bumped body. 


(b) Ibiden'hed bo^. 


Figure 15.- Variation of with several values of length 

z^tio ^2/^1 location 







